Exact Analysis of Entanglement in Gapped Quantum Spin Chains 
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We investigate the entanglement properties of the valence-bond-solid states with generic integer- 
spin S. Using the Schwinger boson representation of the valence-bond-solid states , the entanglement 
entropy, the von Neumann entropy of a subsystem, is obtained exactly and its relationship with the 
usual correlation function is clarified. The saturation value of the entanglement entropy, 2 log 2 (S+l), 
is derived explicitly and is interpreted in terms of the edge-state picture. The validity of our 
analytical results and the edge-state picture is numerically confirmed. We also propose a novel 
application of the edge state as a qubit for quantum computation. 
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Entanglement properties of quantum spin systems 
have been attracting much attention in quantum infor- 
mation theory and condensed matter physics. The entan- 
glement entropy (EE), the von Neumann entropy of the 
reduced density matrix of a subsystem, is a measure to 
quantify how much entangled a many-body ground state 
is. Recently the EE has been used to investigate the na- 
ture of quantum ground states such as the quantum phase 
transition and topological/quantum order [3, HI, 01, 13, [H . 
Vidal et aL[l| conjectured that the EE of a large block of 
spins in gapped spin chains reaches saturation while that 
in critical spin chains shows a logarithmic divergence. 

In this Letter we study the EE of gapped quantum spin 
chains with arbitrary integer-spin. After the Haldane 
conjecture that integer-spin antiferromagnetic Heisen- 
berg chains have a finite gap0, 0], Affleck, Kennedy, 
Lieb and Tasaki (AKLT) proposed the valence bond solid 
(VBS) state which enables us to understandground state 
properties of the Haldane gap systems d, [i| . The VBS 
is now attracting renewed interest from the viewpoint 
of quantum information theory. For example, universal 
quantum computation based on the VBS states has been 
proposed [To( | . 

While the entanglement properties in £> = 1 VBS has 
been extensively studied in Il3j |, we investigate the 
EE in generic VBS states with arbitrary integer spin 
S. We stress that there exist not only S = 1 antiferro- 
magnetic Heisenberg chains [HI, [H| but an S = 2 chain 
(MnCi3(bipy))[l6| in which the presence of the Haldane 
gap has been experimentally confirmed. We give the ex- 
act form of the EE in generic VBS states in this Letter. 
Then we explicitly confirm that the part of the conjecture 
proposed by Vidal et al. is true for all integer-spin VBS 
chains. The relationship between the EE and the corre- 
lation function is clarified and the physical meaning of 
the EE in gapped models is established. We also make a 
comparison between the analytical results for VBS chains 
and the numerical results for higher-spin antiferromag- 
netic Heisenberg chains. The obtained results indicate 
that the edge-state picture is valid not only for 5 = 1 



Haldane chains but also for all the other integer spin-S* 
chains. This is a typical consequence of the non-trivial 
topological and/or quantum orders, where characteristic 
features are hidden in the bulk and appear only near the 
boundaries and impurities [17| ■ We also discuss a poten- 
tial application of the edge states as qubits for quantum 
computation. 

Let us start with the Schwinger boson representation 
of generic VBS states. The spin operators are repre- 
sented by the Schwinger bosons as — a^bj, Sj~ = bjdj, 

ij and bj satisfy 
, ,. , j tJ with the all the other commu- 
tators vanishing [la ]. To reproduce the dimension of 
the spin-5 Hilbert space at each site, we must impose 
the constraint that the total boson occupation number 
OjGj + bjbj = 2S. Using the Schwinger boson represen- 
tation, the spin-S 1 VBS state with two spin-5/2's on the 
boundary is written as 



and 



Sf = 



(ajdj — bjbj)/2, where 
[biM] = Si. 



|VBS) = n>t6t +1 - 6ta) +1 ) s |vac), 



(1) 



where j = 
end sites. 



= 1,2, 
Bi. 



, N are bulk sites and and N + 1 are 
i\b\ — b\a\ is a creation operator for the 



valence bond between i and j [19( . The VBS state ([Tj) is 
a zero-energy ground state of the following Hamiltonian: 



H 



N-l 2S 



E 



1 J=S+1 



ArP, 



where the projection operator P/j + i projects the bond 

spin Jj. j+i = Sj+Sj+i onto the subspace of magnitude J. 
Here the coefficient Aj can be an arbitrary positive value. 
The boundary terms describing interaction between spin 
S/2 and spin S are explicitly written as 

3S/2 3S/2 

tto.i = ^ B,jPq 1: ttn.n+i = ^2 b jPn,n+\i 

J=S/2+l J=S/2+l 
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with Bj > 0. In order to calculate reduced density ma- 
trices, it is convenient to introduce a spin coherent state. 
For a point Cl = (sin6*cos</>, sin^sin^, cos#) on the unit 
sphere, the spin coherent state at each site is defined as 

where = (cos(6»/2)e^/ 2 , sin(6»/2) e -^/ 2 ) are spinor 

coodinates. Here we have already fixed the U(l) gauge 
degree of freedom since it has no physical content. Using 
the trace of any operator O is written as TrC = 
2f±l/dft<ft|0|ft>. " 

Let us now calculate the EE of a block of L contiguous 



bulk spins in the VBS state (p}. For the density ma- 
trix of our ground state p = |VBS)(VBS|/(VBS|VBS), 
the reduced density matrix of the block of L contigu- 
ous bulk spins is defined as pl — r Tr-g L p. Here Bl is 
a block of L spins and Bl is its complement. The EE 
Sl = —TrB L pLlog 2 pL is determined by eigenvalues of 
PL- Suppose that the block of L contiguous spins start- 
ing from site k and stretching up to k + L — 1, where 
k > \ and k + L - I < N (Fig. 1). To obtain the re- 
duced density matrix pl, we take the trace over the sites 
j = 0, 1, k - 1 and j = k + L, ...,N, N + 1. Using the 
spin coherent state representation, pl is formally written 
as 



KUjeBr. Tito C-^^fUl^L (^%^) S QlP^lVBS^VBS^IQ.P^ 



where boundary operators and a block of VBS state 
with length L are defined as Qfe = (itfc— i&fc — 
w fe _ia fc ) 5 , Pk+L-i = (a-k+L-iVk+L - b k+ L-iu k+L ) S and 
|VBS L ) = nSfc" 2 ^^! ~ 6]a] +1 ) S |vac L ), respec- 
tively. Here we have already used the following relation: 
(0\a s - l b s+l \n) = yf(2Sy.u s - l v s+l . In Eq. ©, the inte- 
grals over Qf.-1-i {I = 1,2, — 1) can be performed 
by regarding flk~i as a polar axis. The same holds for 
&k+L+m ( m — 1)2, —,N ~ L — k + 1). After integrat- 
ing over these variables, we immediately notice that the 
reduced density matrix pl does not depend on both the 
starting site k and the total length N. The same prop- 
erty for S = 1 VBS has been proved in [l3[ by another 
approach, i.e. using the special property of maximally 
entangled states. The coherent state approach, however, 
allows us to generalize this result for more complicated 
cases. For example, we can also prove that the EE does 
not depend on the whole size of a VBS state on a two- 
dimensional Cayley tree 20] by using the coherent state 
representation. 

Since the reduced density matrix does not depend on 
both k and N, we can set N = L without loss of gener- 
ality. The following remarkable property makes it easier 
to calculate the EE of L contiguous spins: Sl = S L = 
—Tr-g L Pi\og 2 pi, where p L = Trg^p. One can easily show 
this by using the Schmidt decomposition. Then all we 
have to do is to obtain the eigenvalues of the reduced 
density matrix of two end spin-S'/2's p L = 

/jtiflg ( 1 - fifc 2 " fc+1 ) S P t Ql +1 |v aC )(va C |PoQ £+1 

(3) 



where | vac) (vac | = |0) (0| O |0) i+ i(0|. The state P T |0) 
in the numerator of ([3]) is explicitly given by (a vl — 
b o ul) s \0). From the definition of the spinor coordinates, 
we notice that (u,v) changes to (iv*,— iu*) when we 
change variables from (9, 4>) to (tt — 9, <fi + tt). Then we 
can rewrite -Pj|0)o as (— i) s y/Sl | — £li)o. In the same 
way, Q l +i\0)l+i can be rewritten as ^^/Sl | — 
Substituting these results into Eq. |3]) and changing the 
variables of integration from £lj to — f2j (j = 1, 2, ...,£), 
we obtain p^ = 

r(T-f L \ t-tL—1 j i-n fc -n fc+ i \S 

J Ulj=l 4tt ) llfc=l I 2 ) 

(4) 

Now the physical meaning of p L is quite clear. Eq. (j4|) 
can be regarded as a correlation function between density 
matrices |fii)o(^i| and \&l} l+i(&l\- More precisely, the 
matrix elements of p L are completely determined by the 
two-point correlation functions of the corresponding one- 
dimensional classical statistical model [19|. This can be 
checked by using the binomial expansion of Pq and Ql+i- 
While this interpretation enables us to understand the 
relation between the EE and the correlation functions, it 
is more convenient to use the form ((4]) for the calculation 
of the EE. 

From now on, we follow Ref. [ll[ and obtain the eigen- 
values ofp L . In Eq. gj), T feife+ i = ( J_g*L%±i )S actg ag a 
transfer matrix of the corresponding classical statistical 
model. Expanding Tk.k+i in terms of Legendre polyno- 
mials and using the addition theorem for spherical har- 
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monies, the transfer matrix can be rewritten as 

s i 



-in 



k+l 



s- 



j5>(0 E Yr(n k )Yr(Ci 

1=0 rn=-l 



k+l) 



(5) 



with A(Z) = (-1) ; 5!(5 + l)!/[(5 - Z)!(5 + I + !)•]• Then 
we substitute ([5]) into (j4|), recall the orthonormality of 

spherical harmonics, i.e., / ciOy, m (0)^7 i '(0) = 8 w 5 mm ' ', 
and obtain 



4tt 



PL = 



(5 + 1) 



1=0 



m=—l 



where irreducible Z-th order spherical tensor operators 



m(m) /_ 
2-S/2+1 



-Z, — Z + 1,...,Z) are defined as Tf 71 ' 



J dn|f2>yj rn (n)(n|. We should note here that 
T^™^ acts on the Hilbert space of the left-end spin-5/2 
while (Tj)^ acts on that of the right-end spin-5/2. Let 



us now introduce the following formula found in [11 

l 



E pi° 



(m) 



(6) 



m=—l 



where So and 5l+i denote the left-end and right-end 
spm-5/2's, respectively. Here Ij(X) is a j— th order poly- 
nomial in X and determined by the following recursion 
relation 



2j + 3 



j 2j + 3,S-j+U* 



j + l 2j-lV5 + j 

with I (X) = Ji(X) = ^ (5/2+1)3 • Tnc isotropic two 

site tensor operators Ij (So -Sl+i) (j = 0, 1, S) are mu- 
tually orthogonal with respect to the trace inner product 
Tio, L+i(Ij h)- Since Eq. (J6j> is completely determined 
by the polynomials in So ■ <Sx+i, the reduced density ma- 
trix p£ is diagonal in the basis which diagonalizes the 
total spin operator Jo,l+i = So + Sl+i- Therefore, the 
eigenvalues of are given by 



pl(J) 



SfS , 
2l2 +1 



where J(= 0, 1, 2, S) is a magnitude of the total spin 
and each Pi(J) is (2J + l)-fold degenerate. Finally, the 
EE of a block of £ contiguous bulk spins is explicitly 
written as 

s 

5 L = -2(2J + l) Pi (J)log a p £ (J). 

7=0 

Since the reduced density matrix approaches a (S + 
l) 2 -dimensional identity matrix in the thermodynamic 



g 3 



□ □□□□□□□□ 



A AAAAAAAAAAA 




5co(2) 



2 4 6 8 10 12 14 

Number of Sites - L - 

FIG. 1: Fhe EE for S = 1(«), 8 = 2(o), S = 3(a), S = 4(a), 
S — 6(B), and S — 8(D) VBS chains as a function of the 
block size L. The broken line indicates the saturation value 
500(2). 



limit L — ► co, we can see that Sl < 21og 2 (5+l) = Soo(S) 
and approaches this upper bound exponentially fast in 
L. This saturation can be observed in Fig. [T] where 
the EE Sl for various spin-S 1 VBS chains are plotted as 
a function of the block size L. Here we confirm that 
the conjecture proposed by Vidal et al. is valid for all 
integer-spin VBS chains. 

Next we make a comparison between the above re- 
sults for the VBS chains and numerical results for the 
integer-spin Heisenberg models. Since S = 1 systems 



have recently been extensively studied 24( , we study nu- 
merically the EE and the energy spectra of the 5 = 2 
Heisenberg model and its continuous deformations. One 
of the simplest 5 = 2 Hamiltonian which interpolates 
between these two models can be written as H — 



N 



z j Si ■ S l+ i + a j -(Si ■ S i+ i) 2 + gg(<% ' Si+i) 3 + — 



where a = and a = 1 correspond to the Heisenberg 
model and the 5 = 2 AKLT model, respectively. 

The edge-state picture in general 5 Haldane systems 
allows us to interpret the spectra as follows. The low- 
lying (5+1) multiplets have (2S to tai + l)-fold degener- 
acy in each sector when the system has open boundaries. 
These generalized Kennedy triplet states are almost de- 
generate, and are completely degenerate at the AKLT 
point. This can be understood from the VBS picture. 
It would be valid for the Heisenberg model by some re- 



sults from numerical calculations [22], |23j. Let us now 
show that the ground state properties remain unchanged 
through the adiabatic continuation from the AKLT to the 
Heisenberg model. Fig. [5] (a) shows the a dependence 
of the energy gaps between the ground and the lowest 
two excited states computed by exact diagonalizations of 
the system of N — 10 sites with periodic boundary condi- 
tions. There is no level crossing between the ground state 
and the first excited state, which suggests that the low- 
energy behaviors of the system are adiabatically equiva- 
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FIG. 2: (a)Energy gaps between the ground and the lowest 
two excited states in the system of N = 10 sites with periodic 
boundary conditions. (b)The EE of the S = 2 periodic N = 
6, 8, 10 Heisenberg model and its continuous deformations. 



lent with each other in this parameter region. 

Finally, let us discuss the EE in our system. The ob- 
tained results of the EE from exact diagonalizations are 
shown in Fig. [2] (b). The EE at the AKLT point has 
a tendency to converge to the value <Soo(2) = 21og 2 3 = 
3.16993 as the system size increases. This value coincides 
with our analytically calculated one with open boundary 
conditions (See Fig. [1]). The lower bound of the EE in the 
calculated region is given by S > 2 log 2 3, and this is the 
contribution from the boundaries of the system created 
by taking partial trace over the subsystem. This lower 
bound is equal to the EE at the AKLT point. Taking 
the edge-state picture into account, we can see that this 
lower bound is closely related to the number of degrees of 
freedom emerging at the edge. In other words, if we can 
prepare a sufficiently long spin-5 1 VBS chain with open 
boundaries, each edge state behaves as a free spin-S/2. 
This (S + l)-level system can be used as a qubit(qudit) 
for quantum computation by locally applying a magnetic 
field at the edge. The EE provides a typical measure for 
the quantum resources. We should note here that the 
EE has contributions not only from the edge state but 
also from the bulk except for the AKLT point. In this 
meaning, the AKLT point is a special point since the EE 
has a contribution only from the edges created by tak- 
ing partial trace. This fact is related to the minimum 
correlation length at the AKLT point. It is also inter- 
esting that the EE at the S = 2 AKLT point takes the 
minimum value. A similar behavior has been observed in 
the case of S = 1 [IS]- Thus, we can conjecture that the 
EE takes a minimum value at the AKLT point in general 
SU(2)-invariant models with integer-spin S as far as the 
edge-state picture is valid. 
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